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Abstract. Narrow-band radio spikes have been recorded during a solar flare with unprecedented resolu- 
tion. This unique example allows us to study the effect of low resolution in previously published peak-flux 
distributions of radio spikes. We give a general, analytical expression for how an actual peak-flux distri- 
bution is changed in shape if the peaks are determined with low temporal and/or frequency resolution. It 
turns out that, generally, low resolution tends to cause an exponential behavior at large flux values if the 
actual distribution is of a power-law shape. The distribution may be severely altered if the burst-duration 
depends on the peak-flux. The derived expression is applicable also to peak-flux distributions derived 
at other wavelengths (e.g. soft and hard X-rays, EUV). We show that for the analyzed spike-event the 
resolution was sufficient for a reliable peak flux distribution. It can be fitted by generalized power-laws 
or by an exponential. 

Key words. Acceleration of particles — Methods: statistical — Sun: flares — Sun: corona — Sun: radio 
radiation 



1. Introduction 

Statistical flare models envisage the flare process as an 
ensemble of sub-processes and do not focus on the sin- 
gle constituent processes. The most prominent of these 
global models are the so-called Cellular Automaton mod- 
els (Lu & Hamilton 1991; Lu et al. 1993; Vlahos et al. 1995; 
Georgoulis & Vlahos 1996; Galsgaard 1996; Georgoulis 
& Vlahos 1998; MacPherson & MacKinnon 1999; Isliker 
et al. 2000, 2001). They assume flares to be fragmented 
and stochastic processes, and so the need for a compari- 
son through statistical quantities is present. Furthermore, 
there is direct evidence that flares are fragmented to 
some unknown level (deJager & deJonge 1978, Benz 1985, 
Aschwanden et al. 1990), and that flares really are stochas- 
tic processes (Isliker & Benz 1994; Isliker 1996; Ryabov et 
al. 1997, Veronig et al. 2000). 

Observed peak-flux distributions of flares and flare 
fragments are used to test statistical flare models. Systems 
in a state of self-organized criticality (SOC) lack any char- 
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acteristic scale and thus show power-law distributions. 
Stochastic growth in uncorrelated regions of instability 
yield log-normal peak-flux distributions found in inter- 
planetary type III bursts (Cairns & Robinson 1997). Open 
driven plasmas yield burst-like pulses distributed expo- 
nentially at large fluxes (Robinson et al. 1996). 

There exists a number of observational studies of peak- 
flux distributions of flare-related (non-thermal) emissions, 
in the hard X-ray range as well as in the radio range 
(see references in Aschwanden et al. 1998). In this arti- 
cle, we will concentrate on narrow-band, millisecond spike 
events in the radio range. Robinson et al. (1996) have de- 
termined a peak-flux distribution from a single-frequency 
spike observation, which they find to be exponential for 
high-flux values. Aschwanden et al. (1998) have analyzed 
some spike events with poor temporal resolution (typically 
one measurement point per spike event), and found expo- 
nential distributions. Meszarosova et al. (2000) analyzed 
single-frequency measurements of spikes. They found ex- 
ponential and power-law distributions, the latter being 
very small in extent, however (much less than one decade). 
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Such analyses raise a question: since every measurement 
is only at discrete points in time and frequency, a de- 
tected peak in an observation is in general not identi- 
cal with the true peak which would be seen if continu- 
ous recording were feasible. The detected peak (termed 
pseudo-peak in the following) is likely to be further away 
from the true peak, the lower the time- and frequency- 
resolutions. Therefore the derived peak-flux distribution 
must be expected to be biased. Furthermore, peak-fluxes 
in the radio-range have mostly been determined at fixed 
frequencies, neglecting completely the fact that one may 
be far from the true peak in the frequency-direction. The 
question therefore is what bias in a peak-flux distribution 
must be expected due to finite and possibly low resolution 
in time and/or frequency. 

We give an analytical expression of how a given true 
peak-flux distribution is changed when determined with 
finite time- and frequency- resolution (Sec. 2). The expres- 
sion is independent of the wave-length range under con- 
sideration, it can also be applied e.g. to soft and hard 
X-rays, or to EUV. We then present the peak-flux distri- 
bution of narrow-band radio spikes in an event measured 
with unprecedented high time and frequency resolution 
(Sec. 3). We will discuss the peak-flux distribution of the 
spike event with the introduced statistical theory, as well 
as the distributions reported in the literature, which are 
subject to poor resolution, or even derived without fre- 
quency information. We will address the question of how 
good the time and frequency resolution must be (in terms 
of duration and bandwidth of the events), in order that 
the detected peak-flux distributions are near to the true 
ones, and what is to be expected if the frequency infor- 
mation is not available or not taken into account (Sec. 
4). 

2. The biasing of peak-flux distributions through 
finite resolution 

In this section, we will establish the analytical expres- 
sion which relates the peak-flux distribution of the true 
peaks to the distribution of the observed ones with given 
time and frequency resolutions, which we will refer to as 
pseudo-peaks. We start by making the following defini- 
tions: 

Let a be the amplitude of a pulse (e.g. spike), and assume 
a pulse-shape of the form 

$(v,t,a) = a ■ §„(vg{a)) ■ § t (th(a)) = a - s (1) 

with $i/(0) = 3>t(0) = 1, and g(a) and h(a) functions of 
the amplitude a (for instance g(a) = 1, a, or 1/a), and 
where we have introduced the abbreviation s — s(v,t;a). 
Such a pulse-shape is reasonably general, but of course 
not completely so: we assume that the profiles in the fre- 
quency and in the time direction are independent, and par- 
ticularly, it is assumed that the amplitude causes merely 
a scaling of the duration and/or bandwidth, and that 
there is no dependence on an additional, hidden, possi- 
bly stochastic parameter. 



In statistical language, we consider the ideal measure- 
ment of the true peak-fluxes as the outcome of a random 
variable A, with probability density pa{o) (oi < o < 0-2) 
(which, in other words, is the normalized peak-flux distri- 
bution). Analogously, the measurement of pseudo peak- 
fluxes (the ones subject to finite temporal and spatial res- 
olution) is considered as the outcome of a random variable 
R, with probability density pn(r). The question is what 
the relation between pa{o) and Pn{r) is, i.e. we need to 
find the connection between A and R: 

Measuring the pseudo-peaks can be viewed as choos- 
ing a random point (N, T) in the ^-t-plane within a cer- 
tain rectangle around the true peak, whose side-length 
are tj and r v (the temporal and spatial resolutions, re- 
spectively), and reading out the flux- value at this point. 
Hence, the uniform probability for choosing a random 
point (N, T) in the rectangle of the f-i-plane is trans- 
formed into a probability distribution on the flux axis: for 
a given pulse with amplitude A = a, the random point 
(N, T) transforms into a random point a ■ S on the flux 
axis through the pulse-shape, a ■ S = a ■ $>v(N) $t(T), 
which will give the probability distribution ps(s) of the 
random variable S. More generally, we assume the pulse- 
shape to depend on the amplitude of the pulse (Eq. 1), so 
that S is given as 

S = * v (N-g(a))-$ t (T-h(a)). (2) 

Its probability distribution is conditionally dependent on 
the amplitude, ps(s\a), and its range is Si(a) < S < 
S2(a) = 1 (ps(s\a)ds denotes the conditional probability, 
i.e. the probability for S to assume values in [s,s + ds], 
given that A is known to assume the value a). The random 
point on the flux-axis a ■ S is the pseudo-peak flux, so that 
the relation between the true and the pseudo peak-flux is, 
in terms of random variables, 

R = A ■ S (3) 

and the wanted pseudo-peak flux distribution pn(r) is 
given as 

Pn{r) =PA-s{r)- (4) 

Of course, it would be of interest to invert the prob- 
lem, i.e. to derive from given pn(r) and ps(s) the true 
distribution pa{o)- However, this is not possible, as will 
be shown in Sec. 2.4, since pn{r) andps(s) do not contain 
enough information to uncover pa{o)- 

Three tasks are to be worked out now: 

1. The rectangle in the ^-i-plane, out of which a random 
point is chosen, has to be determined (to find the prob- 
ability density of the point (N, T) , which is needed to 
derive ps through Eq. 2). 

2. The probability ps(s | a) has to be derived (through Eq. 
2), since it is an input to Eq. (4). 

3. The probability distribution of the product of the ran- 
dom variables A and S has to be found to evaluate Eq. 
(4), which yields the wanted pr. 



H. Isliker and A.O. Benz: On the reliability of peak-flux distributions, with an application to solar flares 



3 



2.1. The rectangle around the true peak in which the 
pseudo-peak is located 

The probability density pN,r{v,t) for the random point 
(N, T) in the i/-£-plane, at which the pseudo-peak flux is 
measured, is uniform in some region around the true peak 
(this follows from the complete absence of correlations be- 
tween the measurement and the measured). In order to 
specify px,T{v,t), the shape and size of this region have 
to be determined. 

Without loss of generality, we may assume a true peak 
to occur at to = 0, v$ — 0. Since we assume the burst- 
profile to factorize (Eq. 1), we can expect the measurement 
point (A, T) to lie in a rectangle, 



(N,T) E [L V -,L V+ ] x [L t -,L t+ ] 



(5) 



around the true peak at (to, vq) — (0, 0), and we may de- 
rive the time- and the frequency-interval independently. 
Note that we explicitly treat the case where the pulse- 
shape is asymmetric. If the pulse shape is symmetric, 
then the rectangle is simply given as [— r„/2, r„/2] x 
[— Tt/2, Tt/2], as follows straightforwardly from symmetry 
considerations. 

We start with treating the time-direction. Assume a 
true peak to be located at to = 0. In the measurement 
procedure, a grid of time-points i, is put onto the i-axis 
(with ti — ti—i = Tt, Vi, with tj the time resolution), 
which is randomly positioned relative to to = 0. At one 
of these points, say at ti, the measured flux will be high- 
est, and a pseudo-peak is detected. This means that (i) 
a$ t (U) > a$ t (U + n), and (ii) a$ t (U) > a$ t (U - n) (for 
now we assume the burst profile to be independent of the 
amplitude a) . Obviously, the pseudo-peak occurance-time 
ti lies in the interval 

U G {t'\$t(t') > $(t' + T t ), and $ t (t') > $ t (t' - T t ), 

and \t'\ < Tt } (6) 

The left boundary L t - is given by the equation $ t (L t _) = 
$>t(Lt-+Tt), with allowed values of Lt- in [—Tt, 0], and the 
right boundary Lt+ is analogously given by the equation 
$ t (L t +) = §t{Lt+-n) in the range L t+ G [0, r t ]. In order 
that these solutions are unique, we need to demand that 
the pulse-shape is convex ($" < 0), i.e. strictly increasing 
until the peak, and then strictly decreasing with time, 
which is reasonable for a pulse-shape. If L t -+T t is inserted 
into the equation for L t + , then it is seen that it solves this 
equation, whence 



L t+ = L t - + r t 



(7) 



- the size of the interval out of which a random point 
is chosen in a measurement is Tt, the time- resolution. Of 
practical interest in the following will be that 



MLt-) = MLt+) 



(8) 



So far, we have omitted the scaling factors h(a) and 
(/(a). L t -, for instance, would actually have to be deter- 
mined by the equation § t (Lt-h(a)) = $t((Lt_ — T t )h(a)) 



instead of the one above, and in general one expects 
the solution L t - to depend on a. Under quite general 
assumptions, however, L t - is independent of a, e.g. for 
Gaussian, exponential, or power-law pulse-shapes. All 
these examples are pulse-shapes with the general form 
a$ t _ = a<f>o(c + 4>-(th(a))) for the left branch, and 
a$ t+ = a(j>o(c + <fi+(t h(a))) for the right branch, where </>o 
is any invertible function (e.g. an exponential), c a con- 
stant, and <p- and <p + may be different, but homogeneous 
of the same degree (i.e. <f>i(bt) = b K 4>i(t), for i = — , +, and 
k a constant), e.g. a power-law. 

Completely analogously, L v+ and L v - are determined, 
and again L u+ = L v - +t„ (t v denotes the frequency reso- 
lution) and <& V (L V -) = § V (L V+ ) hold. With the determi- 
nation of the rectangle in which the random point (N, T) 
lies, the probability density pn,t follows immediately as 



PN,r{v,t) = p N {v)PT{t) 



1 

T t T v 



(9) 



since it is uniform, as explained at the beginning of this 
subsection. 



2.2. The determination of Ps(s\a) 

Instead of the probability density ps{s\a), we will de- 
rive the cumulative probability distribution Ps(s\a) := 
// Ps(s'\a) ds' , since in general the density has a singular- 
ity at the peak (s = 1) of the pulse. Ps(s\a) is given implic- 
itly by Eq. (2) as a function of the pulse-shape. Starting 
from the definition of Ps(s\a) and inserting Eq. (2), we 
have 



Ps(s\a) = prob S < s given that A = 

= prob[$„(7V • 9 {a)) ■ $ t (T • h(a)) < s 
= prob (T, N) E 

{(t,v) I $ v (vg(a)) ■ $ t (th(a)) < s, and 
Lv- < v < L v+ , and L t - <t< L t +} 

= J J PT{t)PN{v)dtdv 



(10) 



$„(i/- fl(a)) * t (f-%)) < 3 
L v - <v< L v+ 
L t -<t< L t+ 

where pT(t)pN(v) is given by Eq. (9). 

Because of the possible asymmetry of the pulse-shape, 
it is necessary to treat the four quadrants separately. In 
each quadrant then, the integration limits in Eq. (10) im- 
ply four different cases, depending on the value of s. If 
in the four quadrants we denote the respective inverses of 
the pulse-shape by $t+, &t-> an d 3^+, $v~, and if we 
write L t where one can insert either L t + or L t - without 
changing the numerical values of the respective expres- 
sions (due to Eq. 8), and analogously L v , then, in the first 
quadrant, we have (a substitution t := <& v (v)$t(t) allows 
to calculate the t- integral in Eq. 10): 
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0. for s < $ v (L v g{a))$ t {L t h{a)) 

P s (s\a) = (11) 

1. for <& u {L u g{a)) < s, and $ t (L t h(a)) < s, and s < 1 
\L v+ L t+ \ 



P s (s\a) 



1 1 

g{a) h{a) 



(12) 



II. tov$ v (L v g(a))<s<<S> t (L t h(a)) 
\L v +L t+ \ 



Ps(s | a) 



T v T t 
1 1 



- rt g{a)h{a) 



/ 

*W(* t (I, t \(a))) 



1 $ v (v) 



1 



^ v+ \<S> t {L t h{a)) ) g(a) 



(13) 



III. for $ t (L t h(o)) < s < <& V (L V g(a)) 
\L v+ L t +\ 



Ps{s | a) 



1 1 



- r * g{a)h{a) 



L„ + g(a) 



(14) 



IV. for § v {L v g{a))$ t {L t h{a)) < s, and s < § v {L v g{a)), 
and s < $ t (X t ft-(a)) 



Ps(s|a) = 



\L v+ L t+ \ 



T v T t 
1 1 



™ g(a)h(a) 



L v+ g(a) 



■'+(*t(^t M«0)) 



1 Mf) 



V. for 1< s 



Pg(s|a) = 



$t(L t /i(a)); g{a) 



2.3. The probability distribution of the product A ■ S 

In the last step, we determine the probability distribution 
Pu of the pseudo-peaks, i.e. the probability distribution of 
the product A ■ S (Eqs. 3 and 4). The case of multiply- 
ing two random variables which are independent and have 
infinite range is found in standard textbooks. However, 
since both, A and 5, have finite range, and since S is con- 
ditionally dependent on A, it is worthwhile to give the 
respective formulae. Essentially, the probability distribu- 
tion of the product A ■ S equals pA,s{ a , s ), integrated over 
the subregion of the region [01,02] x [si(a), s 2 (a)] where 
as < r is fulfilled (pA,s( a , s ) is the joint probability distri- 
bution of A and S). Using that pA,s( a , s ) = Pa{o)ps{s\o), 
we explicitly have 

P R (r) = prob[i? <r] = prob[A ■ S < r] 
PA,s{ a , s ) dads 



(15) 



(16) 



The non-trivial cases are for intermediate ranges of s, i.e. 
for $ v (L v g(a))$> t (Lth(a)) < s < 1. The formulae for the 
other three quadrants are gained completely analogously 
to the given ones, by just correspondingly interchanging 
the indices t+, t—, v+, v— (it is for this purpose that we 
had to write absolute values for all the appearing terms). 
If the pulse-shape is symmetric in v and t, then the con- 
tributions of the four quadrants are equal. 



// 



as < r 
ai < a < a 2 
Si(a) < s < S2(a) 



Pa(o) da 



p s (s\a)ds (17) 



a\ < a < a 2 si(a) < s < 52(a) 

a < r/sx(a) s < r/a 

where siai < r < s 2 a 2 . (As mentioned in Sec. 2.2, ps(s\a) 
might have a singularity at s = 1, and we must therefore 
use the cumulative probability distribution Ps(s\a), which 
is achieved if also for R we use the cumulative probabil- 
ity distribution P R {r), as done in Eq. (17), and identify 

J? p s (s>\a)ds> = P s (f3\a) - P s (a\a).) 

In the following, we write again L t if inserting L t + 
or L t - does not change the respective numerical val- 
ues, and correspondingly L v is used instead of L v+ 
or L v -. For intermediate values of r, i.e. for a\ ■ 
<&v{Lvg(ai))<$>t{Lth{ai)) < r < a 2 , four non-trivial cases 
of integration limits turn out to exist, depending on the 
values of r. We have: 

0. for r < ai • ^ v {L v g(a 1 ))<^ t (L t h(a 1 )) 

P R (r) = (18) 

1. for r < a 2 ■ § V {L V g{a 2 )) $t{L t h{a 2 )), and r < ai 

a<r/<S>,4L v g(a))<S> t (L t h(a)) 

P R (r)= J dap A (a)P s (r/a\a) (19) 

ai 

II. for ai < r < a 2 ■ <$> V (L V g(a 2 )) <5> t (L t h{a 2 )) 

a<r/<S> v (L v g(a))<l> t (L t h(a)) 

P R (r)=P A (r)+ J da PA (a)P s (r/a\a) (20) 

r 

III. for a 2 ■ $„(L„ g{a 2 )) <J> t (L t h{a 2 )) <r <a x 

a,? 

P R (r) = J dap A (a) P s (r/a | a) (21) 

ai 
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IV. for a 2 g(a 2 )) $t(L t h(a 2 )) < r, and a\ < r, and 

r < a 2 



a-2 

P R (r) = P A {r) + j dap A (a) P s (r/a \ a) 



V. for a 2 < r 
P R {r) = 1 



(22) 



(23) 



Inserting Ps(s\a) (Eqs. (11) to (16)) into the formulae for 
Pr(t) yields the desired pseudo peak- flux distribution. 

2.4. Inversion is not possible 

To derive the true peak flux distribution from a given 
pseudo-peak flux distribution (the inverse problem), we 
have to proceed as follows: from Eq. (3) we find 



A = R/S, 



(24) 



and we can uncover the true distribution from the pseudo 
one analogously as we had proceeded in Eq. (17): 

P A (a) = prob[A < a] = prob[i?/5 < a] 

= I Y PR,s(r,s)drds 



r/s<a 



J J p R {r\s)p s {s)drds 



(25) 



r/s<a 



so that we need to know the conditional probability 
Pfi( r l s ) 01 R given that S is known. That R is not inde- 
pendent of S is evident from the fact that R is equivalent 
to A ■ S, and one explicitly finds, if is assumed to 

be known, that 



P R (r\s) = P A (r/s) 

What is measured, however, is 



Pn(r) 



PR(r\s)p s (s) ds, 



(26) 



(27) 



si<s<s 2 



the distribution of R irrespective of the value of S, so that 
all the conditional information is lost. In other words, to 
uncover Pa{o) one has to know Pn(r\s), which is essen- 
tially equivalent to knowing the true peak-flux distribu- 
tion Pa{o) (see Eq. 26), and which would be feasible only 
in continuous observations. Whence it follows that uncov- 
ering the true distribution from the measured one is not 
possible. 

3. The peak-flux distribution of solar narrow-band 
spikes 



Table 1. Fits of different functional forms to the peak-flux 
distributions of the event 1982/06/04 (see Sec. 3), where 
the four data-sets are distributions derived from a) peaks 
in 2D plane, r t = 2ms, t„ = 1 MHz (data-set a); b) peaks 
at v = 362 MHz, r t = 2 ms (data-set b) ; c) peaks in 2D 
plane, r t = 100 ms, t„ = 1 MHz (data-set c); d) peaks at 
v = 362 MHz, r t = 100 ms (data-set d). If the \ 2 approved 
a fit, then the power-law index or an 'o.k.' is stated, else 
a ' — ' is noted. 



data 


nr. of 


ax c 


a(x - b) c 


a(x -b) c + d 


ae cx 


set 


peaks 


c = 


c = 


c = 




a 


59 




-15.6±6.8 


-2.7±7.5 


o.k. 


b 


144 






-14.9±7.3 


o.k. 


c 


38 


-1.2±0.2 


-3.5±6.9 


-0.8±7.9 


o.k. 


d 


76 




-16.3±4.8 


-15.8±6.8 


o.k. 



The solar narrow-band millisecond spike event we ana- 
lyzed was observed by the ETH Zurich radio-spectrometer 
on 1982/06/04, 13:38:41 UT (the event is published and 
described in Giidel & Benz, 1990; Csillaghy & Benz, 1993). <S> v {v) = e~M^) 



The resolution is 2 ms in time and 1 MHz in frequency 
(from 361 to 364 MHz), whereas the spikes have a typical 
duration of 73 ms (FWHM; Giidel & Benz, 1990), and a 
typical bandwidth of 7 MHz (FWHM; Csillaghy & Benz, 
1993). Hence the spikes are well resolved in time. Also in 
frequency, the spikes are resolved although the observa- 
tion range (4 MHz) is smaller than the typical bandwidth, 
as only spikes with peaks unambiguously in the range are 
taken into account. The peaks were determined by looking 
for strong enough local maxima above the noise-level, and 
a constant background was subtracted (representing the 
quiet Sun) . The normalized distribution of the respective 
peak-fluxes is shown in Fig. 1 (solid line, with error-bars). 

Different curves were fitted to this distribution, and a 
X 2 -test was performed to check whether the fits are com- 
patible with the data or not. The fitted curves are: a sim- 
ple power-law {ax c ), two generalized forms of power-laws 
(a(x— b) c , and a(x— b) c +d), and an exponential (ae x ). The 
result is presented in Fig. 2 and summarized in Table 1 
(data-set a) : The peak- flux distribution of the event can be 
fitted by the generalized power-laws as well as by the ex- 
ponential, but not by the simple power-law. The indices of 
the generalized power-laws are subject to large errors (de- 
termined by the bootstrap resampling method) due to the 
large error-bars (the number of peaks is relatively small) 
and due to the relative flexibility of the generalized power- 
laws (three resp. four free parameters and 10 data-points). 



4. Discussion of the empirical peak-flux 
distributions 

To apply the statistical theory introduced in Sec. 2 to the 
narrow-band spike event analyzed in Sec. 3, we have to 
make an assumption about the pulse-shape of the indi- 
vidual spikes. According to Giidel & Benz (1990), and 
Csillaghy & Benz (1993), it is reasonable to assume a 
Gaussian pulse-shape, i.e. 



(28) 
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0.01 - 




0.0001 r 



100 1000 
amplitude 

Fig. 1. Peak-flux distributions (probability densities) of 
the narrow-band spike event 1982/06/04 (see Sec. 3): a) 
peaks in 2D plane, r t — 2 ms, t v = 1 MHz (solid line, with 
error- bars); b) peaks at v — 362 MHz, r t = 2ms (dashed 
line); c) peaks in 2D plane, r t = fOOms, t„ = f MHz 
(dotted line); d) peaks at v = 362 MHz, r t = 100 ms (dash- 
dotted line). (The histograms are drawn by connecting 
the midpoints of the bins, and the bin-widths are such 
that each contains the same number of data-points. The 
amplitude is in SFU.) 

and 

$t(t) = e~*(^) 2 (29) 

with b v = 3.1MHz, and b t = 31ms (implying a FWHM 
of 7.3 MHz, and 73 ms, respectively, as reported by Gudel 
and Benz (1990), and Csillaghy and Benz (1993) for the 
given event). Furthermore, we have to assume a distribu- 
tion of the true peak-fluxes A: we let Pa{o) be a straight 
power-law 

Pa{o) = Ca~ a , ai < a < a 2 (30) 

with a 2 = 1000 (from Fig. 1). Other true distributions can 
be expected to produce qualitatively the same effects as 
reported below for the case of this power-law. 

The crucial parameters are the time-resolution r t and 
the frequency- resolution t„, they determine how near the 
pseudo-peaks are to the true-peaks, on average. Whence, 
in the following parametric study (Figs. 3 to 7), we al- 
ways show the true peak-flux distribution pa(cl), together 
with the pseudo peak-flux distributions for four cases of 
time- and frequency-resolution: (i) good resolution in time 
and frequency, (ii) good resolution in time and a bad one 
in frequency, (iii) bad resolution in time and a good one 
in frequency, and (iv) bad resolution in both time and 



0.01 r 




0.0001 - 



100 1000 
amplitude 

Fig. 2. Fits to the peak-flux distribution (probability den- 
sity) of the narrow-band spike event 1982/06/04 (see Sec. 
3, and also Table 1). The peaks are determined in the 
2D plane, with r f = 2ms and t„ = 1 MHz. Plotted are 
the original data (solid line), and the fits ax c (dashed), 
a(x — b) c (dotted), a(x — b) c + d (dash-dotted), ae x (wide 
dash-dotted). (The histogram is drawn and generated as 
described in Fig. 1, and the amplitude units are again 
SFU.) 

frequency (by 'good' we mean r << FWHM, and 'bad' 
means r = 2 FWHM). The cases (ii) and (iv) represent 
also the case of single-frequency observations. 

First, we investigate the case of an amplitude- 
independent pulse-shape, i.e. g(a) = h(a) = 1 (see Eq. 1), 
and we set oi — 1 , a = 2. Fig. 3 shows that for good reso- 
lution both in frequency and time the true and the pseudo 
peak-flux distributions practically coincide. If one or both 
resolutions are low, then the pseudo-peak flux distribution 
is generally near the true one, with a faster fall-off at high 
flux-values, however, i.e. a turning from power-law to ex- 
ponential behaviour. At pseudo peak-flux values smaller 
than ai , a completely artificial, relatively flat extension of 
the distribution appears. Turn-overs at small fluxes in dis- 
tributions detected with low resolution(s) might thus be 
just the effect of an intrinsic low-amplitude cut-off in the 
true fluxes. Only high resolution analysis can tell whether 
or not such a flattening is real or not. The turning-over 
effect increases if the resolution decreases. 

We turn now to the case of amplitude dependent pulse- 
shapes. If we let g(a) = h{a) = a -1 , then we get a devia- 
tion from the power-law behavior at small flux values, but 
nice coincidence for intermediate and high values (Fig. 4) . 
For g(a) = h(a) = a 1 the slope (and shape at high flux- 
values) is drastically changed to a stronger fall-off (Fig. 5), 
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Fig. 3. Amplitude-independent pulse-shape: (0) True 
peak-flux distribution (a = 2, a\ — 1, g(a) — h(a) = 1) 
(solid), together with the pseudo peak-flux distributions 
for (i) r t — 2ms, r„ = 1 MHz (dashes); (ii) r t = 2ms, 
t v = 15 MHz (small dashes); (iii) r t = 100 ms, t v = I MHz 
(dots); (iv) T t = 100ms, t v = 15 MHz (dot-dashed). The 
cases (0) and (i) practically coincide. 
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Fig. 4. Amplitude- dependent pulse-shape I: True peak-flux 
distribution (a = 2, a\ = 1, ,g(a) = /i(a) = a -1 ) (solid), 
together with pseudo peak-flux distributions for r t and t„ 
as described in the caption to Fig. 3. The cases (0) and 
(i) practically coincide. 



Fig. 5. Amplitude- dependent pulse-shape II: True peak- 
flux distribution (a = 2, a\ = 1, g(a) — h(a) = a 1 ) (solid), 
together with pseudo peak-flux distributions for Tt and t„ 
as described in the caption to Fig. 3. 

the observed distribution is different even for good time 
and frequency resolution. 

The general behavior demonstrated so far is rather 
independent of the shape of the original, true peak-flux 
distribution: In Fig. 6, we let again g(a) = h(a) = 1, 
but now the power-law index a = 3. The corresponding 
pseudo peak-flux distributions behave analogously to the 
case a = 2 (Fig. 3). 

It is worthwhile noting that the parameter a\ may en- 
hance the effect of turning to exponential behavior at high 
flux values: in Fig. 7, we let g(a) = h(a) = 1 and a = 2, 
as in Fig. 3, but a\ = 30 (cf. the lower cut-off in Fig. 1), 
and the tendency seen in Fig. 3 is enhanced, now. Thus, 
we conclude that the roll-over at large amplitude is more 
serious for a small range of amplitudes. 

Obviously, it is crucial to know whether the durations 
and bandwidths of the pulses depend on the amplitude: 
Csillaghy and Benz (1993) report that the bandwidth of 
spikes sometimes is correlated with the amplitude, some- 
times it is not, and sometimes it is anti-correlated. There 
is, however, no generally holding strong tendency, so that 
we may assume that the peak-flux distribution reported 
in Sec. 3 is near the true one. 

This is also confirmed by the following: we artificially 
worsened the time and the frequency resolution in the 
data of Sec. 3 and compared the respective histograms: 
The dashed line in Fig. 1 is the histogram of the pseudo- 
peak fluxes resulting if the frequency information is com- 
pletely neglected, i.e. one of the frequencies is selected, 
and the peaks are determined as maxima in time-direction 




Fig. 6. Steeper true distribution: True peak-flux distribu- 
tion (a = 3, a\ = 1, g(a) = h(a) = 1) (solid), together 
with pseudo peak-flux distributions for r t and r„ as de- 
scribed in the caption to Fig. 3. The cases (0) and (i) 
practically coincide. 



Fig. 7. Smaller range of amplitudes: True peak-flux dis- 
tribution (a = 2, ai = 30, g(a) = h(a) = 1) (solid), 
together with pseudo peak-flux distributions for r t and t v 
as described in the caption to Fig. 3. The cases (0) and 
(i) practically coincide. 



only. The dotted line is the distribution for the resam- 
plcd observation, using only every 50th data-point in time, 
which yields a spectrogram with 0.1 sec time resolution, 
and keeping full frequency resolution. Finally, the dashed- 
dottcd line is the distribution for the resampled data, with 
again fixing a frequency and neglecting completely the cor- 
responding information. Obviously, the biasing effects are 
smaller than the statistical errors in the distributions, all 
four distributions coincide within the error-bars, and the 
results are fairly independent of the sampling — only the 
fits seem to show different kinds of behavior (Table 1), 
but, as mentioned, the indices of the generalized power- 
laws are subject to large errors. The relative robustness 
(within the statistical errors) of the distribution on under- 
sampling is in turn a signature for amplitude-independent 
pulse-shapes. 

5. Conclusion 

The statistical theory we introduced allows us to predict 
the deviation of a pseudo-peak flux distribution from the 
true peak-flux distribution if the time and frequency res- 
olutions are known. It turns out that in general there is a 
tendency towards exponential behavior at large flux val- 
ues; the more expressed, the lower the resolutions are, in- 
cluding in particular the case of single-frequency observa- 
tions. Only with high resolutions in both frequency and 
time (compared e.g. to the respective FWHM) are the 
detected distributions reliable in the whole range. The de- 



pendence of the pulse-shape on the amplitude (peak-flux) 
is crucial: If the width of a pulse (duration or bandwidth) 
is proportional to the inverse of the amplitude, then a 
strong deviation from the true peak flux distribution will 
result, the distributions will be steepened and completely 
biased in the whole range, even for high resolution in both 
frequency and time. If the width of the pulse is directly 
proportional to the amplitude, then a bias (flattening) ap- 
pears only in the low amplitude range. All the biasing ef- 
fects get stronger with a smaller extent of a distribution. 

A different possible cause for a strong bias at low am- 
plitudes (appearance of a relatively flat part in the de- 
tected distribution) is a possible intrinsic lower cut-off of 
the true amplitudes. Only high resolution data analysis 
can make sure whether a flattening at low amplitudes is 
real or an artifact. 

The example of a narrow-band spike event we analyzed 
is a unique observation with respect to the high tempo- 
ral and frequency resolutions; the spikes are completely 
resolved in frequency and time. Since moreover the burst- 
width seems to depend at most slightly on the amplitude, 
we conclude on the basis of our analytical study that the 
peak-flux distribution we get is reliable. It can be fitted 
by generalized power-laws or by an exponential, but defi- 
nitely not by a simple power-law. 

From our analysis, it follows that the peak-flux dis- 
tributions of spikes reported by Robinson ct al. (1996), 
Aschwanden et al. (1998), and Meszarosova et al. (2000), 
though observed with poor time-resolution or without any 
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frequency information, are near the true distributions, ex- 
cept for the high- flux part, which must be expected to be 
too steep. 

The exponential distribution we find supports the open 
driven plasma model (Robinson et al. 1996). We did 
though not try to fit a log-normal distribution, which in 
view of the relatively large statistical error in the empirical 
distribution, is likely to be also accepted by the x 2_ test, 
so that we cannot exclude the stochastic growth theory 
(Cairns & Robinson 1997). 

Conclusions on whether the distribution we find is 
compatible or not with SOC models are more difficult 
to draw. SOC models arc models for the primary en- 
ergy release, so far they do not include a mechanism for 
radio- (plasma-) emission. Moreover, the peak-fluxes of 
SOC models, which have been analyzed statistically, are 
defined as the peak-fluxes of entire flares, whereas here and 
in the cited articles on radio observations, the peak-fluxes 
of all the flare-fragments are analyzed. For both reasons, 
a direct comparison of SOC models to radio data is not 
possible without substantial new developments. 

The statistical theory introduced here can be applied 
to the peak-flux distributions of all the different kinds of 
bursts, independent of the wavelength at which they oc- 
cur, as soon as the pulse-shape and the functional depen- 
dency of the pulse-shape on the amplitude are (at least 
approximately) known: to radio-bursts (recent empirical 
studies include Aschwanden et al. (1998; type III, deci- 
metric pulsations), Mercier and Trottet (1997; type I)), to 
soft X-rays, EUV, hard X-rays, etc. (see e.g. the review of 
Crosby et al. 1993; Krucker & Benz 1998). 
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